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^ ; ABSTRACT 

2 ■ We present a study of carbon in N159/N160, an H II region complex in the 

c/3 . low metallicity Large Magellanic Cloud. We have mapped this region, which 

comprises four distinct molecular clouds spanning a wide range of star-formation 
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activity, in four transitions: ^^CO (J = 1 — > 0), ^^CO (J = 2 — 1) and 
(J = 4 ^ 3), and [C I] (^Pi ^^Pq)- Combining these data with existing [C II] 

P3/2 ^^Pi/2) observations provides a complete picture of the predominant 
forms of carbon in the gas phase of the ISM. The new CO (J = 2 — ^ 1) 
data show that the complex is immersed in an envelope of extended, low level 
emission, undetected by previous (J = 1 ^ 0) mapping efforts. The ^^CO 
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(J = 2— >1)/(J = 1— >0) ratio in this envelope is ;^ 3, a value consistent with 
optically thin CO emission. The envelope is also relatively bright in [C I] and 
[C II], and calculations show that it is mostly photodissociated: it appears to 
be translucent (Ay < !)• Neutral carbon emission in the complex unexpectedly 
peaks at the quiescent southern cloud (N159S). In the northern portion of the 
map (the N160 nebula), the H II regions prominent in [C II] correspond to 
holes in the [C I] distribution. Overall we find that, while the I[cii]/Ico ratio is 
enhanced with respect to similar complexes in the Milky Way, the I[ci]/Ico ratio 
appears to be similar or reduced. 

Subject headings: Magellanic Clouds — galaxies: ISM — galaxies: irregular — 
radio lines: ISM — submiUimeter — ISM: clouds 



1. Introduction 

The Magellanic Clouds are especially interesting systems for the study of the interstellar 
medium (ISM). Their low heavy element abundances and low dust-to-gas ratios, coupled 
with their proximity and their active star formation, make them ideal laboratories for 
understanding how metallicity and strong radiation fields infiuence the composition and 
structure of the ISM. Furthermore, in several important aspects they resemble primeval 
galaxies: they are morphologically irregular, metal-poor, and very actively forming massive 
stars. It is, however, worthwhile to point out that the Magellanic Clouds are not pristine 
systems. Their irregular morphologies refer mostly to the young populations (massive stars, 
OB associations, and H II regions), probably influenced by recent encounters with the 
Milky Way: the distribution of old stars and overall mass is not irregular (Cioni, Habing 
& Israel 2000). Nevertheless, because of a relatively low time-averaged star formation rate 
the Magellanic Clouds are considerably less developed that the Milky Way. Therefore their 
study will advance the understanding of the interaction between star formation and the 
ISM in primordial systems. 

The H II regions in the N159/N160 complex were first cataloged by Henize (1956) in a 
survey of emission line nebulae in the Large Magellanic Cloud (Davies, Elliot, & Meaburn 
1976 give them the numbers 271 and 284). These nebulae are two of the brightest H II 
regions in the immediate vicinity of 30 Doradus, located ~ 40' to its south. 

In this region the star formation activity appears to be progressing from the (relatively) 
evolved starburst of 30 Doradus, where most of the enshrouding molecular cloud has been 
dissipated (Johansson et al. 1998), towards the quiescent southern CO arm region were 
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little or no star formation is currently taking place (Cohen et al. 1988; Kutner et al. 1997). 

Located at such a transitional place the N159/N160 complex is one of the best studied 
star-forming regions of the Large Magellanic Cloud in molecular and atomic transitions 
(e.g., Johansson et al. 1994; Israel et al. 1996; Stark et al. 1997a; Pak et al. 1998; 
Johansson et al. 1998; Heikkila, Johansson, & Olofsson 1999). 

The N159/N160 complex features three distinct and spatially well separated regions: 
1) the northern region, chiefly associated with the N160 nebula, where massive star 
formation is well evolved and the parent clouds have been mostly, albeit not completely, 
photodissociated and dissipated. 2) The central region, associated with the N159 nebula, 
which is undergoing strong star formation activity but still wrapped in molecular gas. This 
region includes two giant molecular clouds (GMCs) known as N159-east and N159-west 
(N159E and N159W respectively). And, 3) the southern region, featuring the molecular 
cloud N159-south (N159S). This cloud is actually the beginning of the 30 Doradus CO ridge 
region, a ~ 900 pc (~ 1°) long spur of molecular material extending southward of the 30 
Doradus nebula. 

The 30 Doradus CO ridge is the largest concentration of CO in the LMC (Cohen et al. 
1988). Because the ridge is located on the leading edge of the LMC's movement through 
the hot Galactic halo (Mathewson & Ford 1984), its formation has been attributed to ram 
pressure compression of the Magellanic interstellar medium (de Boer et al. 1998; Kim et 
al. 1998). The entire ridge region (including N159S) is quiescent, with little or no star 
formation activity as evidenced by its faint far- infrared (FIR), Ha (Davies et al. 1976), and 
[C II] emission (Mochizuki et al. 1994). 

Because dust is the dominant source of UV extinction, the low dust-to-gas ratio of 
the metal-poor ISM allows the UV radiation to penetrate more deeply into the molecular 
material thereby causing widespread photodissociation and photoionization. The regions, at 
the surfaces of clouds, where the UV radiation dominates the heating and chemistry of the 
ISM are known as photon-dominated or photodissociation regions (PDRs). In these regions 
the UV radiation dissociates most molecular species, but leaves H2 relatively intact because 
of its strong self-shielding (Abgrall et al. 1992; Pak et al. 1998). Therefore in metal-poor 
systems an increasingly large fraction of the molecular gas is expected to be associated with 
C° and C"*", rather than CO. Observations of CO together with its photodissociation and 
photoionization products (C° and C"*") are the key to understand the interplay between the 
radiation field and the chemical and physical structure of molecular clouds in low metallicity 
systems (e.g., Bolatto, Jackson, & Ingalls 1999; Bolatto et al. 2000). This paper presents 
some of these observations. 

What is the dust-to-gas ratio in the LMC and how does it affect its PDRs? Several 
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arguments suggest that the dust-to-gas ratio of a galaxy is proportional to its metallicity 
Z (Franco & Cox 1986), but recent studies find a dependence (Lisenfeld & Ferrara 
1998). Independently of this, the dust-to-gas ratios measured in the Magellanic Clouds 
are ~ 1/4 and ~ 1/8 of the Galactic value (LMC and SMC respectively, Koornneef 1982; 
Bouchet et al. 1985), while their metallicities (measured through the oxygen abundance) 
are Zlmc ~ ^Orion/2 and 2'gMc ~ -Z'orion/5. Because of this relative lack of dust, to achieve 
the same visual extinction Ay in the LMC requires ~ 4 times larger gas columns than in 
the Milky Way. Therefore, PDRs are expected to be on average ~ 4 times more extended 
in the LMC, while the CO cores will be 4 times smaller (e.g., Israel 1997; Pak et al. 1998). 

Observations suggest that CO cores for clouds of similar virial mass are indeed smaller 
in the LMC than in the Milky Way, and most CO is photodissociated. The CO ( J = 1 — > 0) 
emission from LMC clouds is several times fainter than in the Milky Way (Cohen et al. 
1988; Israel et al. 1993; Rubio, Lequeux, & Boulanger 1993; Kutner et al. 1997). Moreover, 
the I[cii]/Ico ratio in the LMC (which measures the fraction of carbon in PDRs) appears to 
be considerably enhanced with respect to the Milky Way (Mochizuki et al. 1994; Poglitsch 
et al. 1995; Israel et al. 1996). 

2. Observations 

In this study we discuss four new datasets on the N159/N160 region: ^^CO (J = 1 ^ 0) 
at z/ ~ 110.2020 GHz, ^^CO (J = 2 ^ 1) at z/ 230.5439 GHz (Fig. [l|), ^^CO (J = 4 ^ 3) 
at z/ ~ 461.0878 GHz (Fig. |), and [C I] (^P^ ^^p^) at z/ ~ 492.1607 GHz (Fig. 0). The 
first two were acquired using the Swedish-ESO Submillimeter Telescope (SEST) at La Silla, 
Chile, during May 1999. The latter two were made with the 1.7 m Antarctic Submillimeter 
Telescope and Remote Observatory (AST/RO, Stark et al. 1997b) at the Amundsen-Scott 
South Pole base, during the austral winter of 1998. Assuming an LMC distance of 52 
kpc [D.M. ~ 18.58; Panagia, Gilmozzi, & Kirshner 1997), the "plate scale" for these 
observations is ~ 15 pc per arcminute. 

To complement these millimeter and sub millimeter- wave datasets, we also use 21 
cm H I and radio continuum data from the ATCA (Australia Telescope Compact Array) 
aperture synthesis survey of the LMC (Kim et al. 1998). Approximately three pointing 
centers were included in this map of the N159/N160 region. The zero-spacing fiux, measured 
using the Parkes telescope, has been added to the ATCA data to obtain reliable column 
density estimates. 

The ^^CO ( J = 1 — i> 0) transition is excited in almost all conditions found in molecular 
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clouds (gas kinetic temperature Tfcj„ ^ 5 K and critical density ricr ~ 1-7 x 10"^ cm~'^) and 
is generally used as a tracer of gas column density, since it is optically thinner than its 
isotopomer ^^CO. The ^^CO ( J = 2 ^ 1) transition requires slightly warmer and denser gas 
to be excited {Tkin ^ H K, ricr ~ 1.1 x 10^ cm~^), but these conditions prevail in all except 
the coldest and most diffuse molecular clouds. The excitation of the ^^CO [J = 4 3) 
transition requires much warmer and denser gas (Tfcj„ ^ 55 K, ricr ~ 10^ cm"'^). Thus ^^CO 
[J = 4 —>■ 3) is excited only in the dense clumps and in the presence of heating sources. 
Finally, the [C I] ('^Pi — ^'^Pq) transition is excited in very general conditions {Tkm ^ 23 K, 
~ 470 cm~^). Because C° is understood to be chiefly produced by the photodissociation 
of CO, this transition is thought to trace the warm PDR that forms in the surfaces of 
molecular clumps exposed to UV (we will, however, explore other alternatives later). The 
precise values of the critical density in all the above cases depend on the actual magnitude 
of radiative trapping effects, the ratio of ortho- and para-H2, and the kinetic temperature 
of the coUisional partners, and are thus only approximate. 

2.1. Millimeter Data 

The bulk of the SEST data has system temperatures in the range 220 — 350 K at 230 
GHz and Tgys ~190 — 250 K at 110 GHz. The maps were acquired during 1999 May, using 
the facility IRAM 115/230 GHz dual SIS receiver. The two frequencies were fed into the 
split high resolution acousto-optical spectrometer (AOS) backend. The resulting bandpasses 
(~ 50 km s~^ for 230 GHz and ~ 100 km s~^ for 110 GHz) are more than adequate for the 
relatively narrow lines and small velocity variation inside N159/N160. The spectra were 
acquired in double beam switching mode, chopping 11.5' in azimuth. A few of the spectra 
showed contamination by emission in the off position, and these positions were reobserved at 
a different parallactic angle. At these frequencies the telescope has HPBW(230)~ 23" and 
HPBW(110)Ri 48", and the main beam efficiency is r/„b(230) ^ 0.50 and r/„b(110) ^ 0.70. 

The map was initially observed on a 24" grid (i.e., approximately full-beam spacing 
at 230 GHz) and the brightest emitting regions were afterwards filled in on a 12" grid 
(half-beam spacing). Two regions selected from the portion of the map showing low 
level emission were reobserved several times simultaneously in ^^CO (J = 1 — > 0) and 
(J = 2 ^ 1) (boxes in Fig. |I|). Pointing and focus on a nearby SiO maser were performed 
before each observing session. The SEST pointing accuracy is 3" rms. 
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2.2. Submillimeter Data 

The ^^CO (J = 4 — > 3) AST/RO observations were made with a system temperature 
Tgys ~ 3500 K. They were acquired on 15 May 1998 using the lower frequency side of 
AST/RO's 460/810 GHz dual SIS waveguide receiver. The [C I] (^Pi ^^p^) observations 
had a system temperature Tsys ~ 1700 K. The region was observed twice, on 1 July and 10 
August of 1998, using the AST/RO SIS quasi-optical receiver. The backend in both cases 
was the 2048 channel low resolution AOS providing a spectral resolution Af ~ 0.4 km s~^ 
over a bandpass of ~ 350 km s~^. The spectra were acquired in position switching mode, 
chopping 15' in azimuth (which is the same as R.A. at the pole). 

The telescope was later determined to have HPBW(461)~ 3.4' ± 0.3 and 
HPBW(492)~ 3.8' ± 0.3, using scans across the limb of the full Moon. The difference 
between the measured beams and the diffraction limited beam size for the telescope was 
caused by a 2° angular misalignment of the tertiary. No successful beam maps were acquired 
during the season, however, and the precise beam shape for either receiver is not actually 
known. Some of the source elongation manifested in the peaks of Fig. ^, for example, may 
not be intrinsic but caused by a distorted beam shape in that particular receiver. 

The forward efficiency for AST/RO was determined from skydips to be 70%, and it is 
assumed to be identical to the main beam efficiency as the telescope is an off-axis, unblocked 
aperture dish. Thus we assume rjmb ~ 0.7 for 461 and 492 GHz. The maps were observed on 
a 30" grid, considerably oversampling the beam. A few of the spectra were contaminated by 
strong noise spikes that occasionally occur in two specific regions of the passband, possibly 
associated with electrical cross-talk between digital and analog signals in the AOS. These 
problematic spectra were discarded. Some grid points are therefore missing, but because 
of the oversampling this does not compromise the quality of the maps. The sensitivity of 
the [C I] map is la ~ 0.1 K km s~^ in the central area, but somewhat worse at the edges 
because there the spatial averaging cannot take advantage of the heavy beam oversampling. 

AST/RO's receivers sit on a static optical bench, located in the warm room below the 
dish, that does not track with the telescope. As a result any optical misalignment translates 
into precession of the receiver beam in the sky, a problem for faint sources that require long 
integrations. AST/RO's radio pointing was carefully controlled during the 1998 season by 
repeatedly observing a set of compact line emitting sources at different elevations and fitting 
a model that compensates for elevation-dependent pointing offsets and beam precession. 
Thus we are confident that the pointing accuracy during our observations was better than 
~ 30" rms. 
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3. Discussion 

3.1. The Molecular ISM as Revealed by CO Observations 

The CO dataset is composed of new ^^CO (J = 1 0), ^^CO (J = 2 ^ 1), and 
^^CO (J = 4 — ^ 3) maps, as well as a preexisting ^^CO (J = 1 ^ 0) map obtained by 
the SEST Magellanic Clouds key program. Figure shows the ^^CO (J = 2 ^ 1) and 
^^CO ( J = 1 ^ 0) SEST data. The four main molecular peaks (Table [^) are apparent in 
both transitions: N160 to the north is a relatively weak and elongated cloud, N159W is 
the strongest peak in the complex, N159E appears to break up into three clumps at this 
resolution, and N159S is a triangular cloud with the smallest (J = 2 — > 1)/(J = 1 — > 0) 
intensity ratio in the complex. The cloud at offset ~ [+2', +5'] is another member of the 
complex, which we will refer to as N160E. This cloud has very narrow lines (FWHM~ 2.6 
km s~^). In addition to the bridge connecting N159S with the northern portion of the 
complex, there is a wealth of extended, faint emission throughout the mapped region. It 
is important to stress this point: even below the 2 K km s~^ contours essentially every 
position in this map shows a significant CO (J = 2 ^ 1) line if neighboring spectra are 
averaged. We will investigate the character of this extended, low level emission in § p.l.2| . 

The overall correspondence of these data with the published ^^CO (J = 1 ^ 0) 
map (Johansson et al. 1994; Johansson et al. 1998) is extremely good, although some 
differences are apparent in the relative intensities and positions of the individual peaks. 
The shift in position is most noticeable for N159S, which at 1.3 mm appears to be ~ 21" 
eastward of its nominal 2.6 mm position. This offset is most probably not real, but caused 
by undersampling. These regions were mapped on a 30" grid at (J = 1 ^ 0), while the 
SEST beam size at 115 GHz is 45" (Johansson et al. 1998). Our own (J = 2 — > 1) data 
are only sampled every 230 GHz beam (24" grid) for N159S (N159E, N159W and N160 
were observed on a 12" grid). The CO (J = 2 ^ 1) peak in N159S falls almost precisely 
halfway between the ( J = 1 ^ 0) samplings, and its position agrees very well with the fully 
sampled ^^CO (J = 1 ^ 0) map. This suggests that the (J = 2 ^ 1) position for this core 
may be better than the ( J = 1 — > 0) location. Excitation gradients within N159S may also 
explain this positional discrepancy. Fortunately, the shifts in the position of the peaks are 
much smaller for the remaining clouds. To maintain consistency with previous work (e.g., 
Heikkila et al. 1999) we use the CO ( J = 1 ^ 0) positions tabulated in Table |I|. Because 
the submillimeter data have much lower angular resolution the precise position of the peaks 
will not greatly affect the ratios in Table |^. 
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3.1.1. Physical Conditions 

What are the physical conditions in the different clouds of the complex? Figure ^ 
compares the ^^CO (J = 1 ^ 0), (J = 2 ^ 1), and the ^^CO (J = 1 ^ 0) observations. All 
three datasets have been spatially smoothed to a common 60" resolution, and the results 
for the regions of interest are summarized in Table |^. The ^^CO/^^CO integrated intensity 
ratio is commonly used to estimate column density. The ^^CO ( J = 2 ^ !)/( J = 1 ^ 0) 
line ratio is sensitive to density in the low density, optically thick limit, and to temperature 
for high densities (e.g., Kaufman et al. 1999). In the optically thin regime this ratio is 
sensitive to temperature, assuming that the emission is thermalized (for an exhaustive 
discussion of CO excitation see, however, Warin, Benayoun, & Viala 1996). 

The centers of the N160 and N159W clouds have the highest ^^GOj^'^CO ( J = 1 ^ 0) 
ratios and therefore the largest opacities, while at N159E this ratio is smaller by a factor 
of ~ 2. Conversely, N160 has the largest ^^CO (J = 2 ^ 1)/(J = 1 ^ 0) ratio and 
therefore appears to be the hottest, while N159S is the coldest cloud in the complex. 
These conclusions agree broadly with a much more detailed multiline excitation analysis 
performed toward the main molecular peaks in the complex by Heikkila et al. (1999), who 
find kinetic temperatures of 20, 25 and 10 K for N159W, N160 and N159S respectively. 

The distribution of the warm and dense gas in the complex is elucidated by Fig. 
Recall that the ^^CO (J = 4 ^ 3) transition requires T ;^ 50 K and n ~ 10^ cm~^ to be 
excited (§0). The angular resolution of the (J = 4 — > 3) observations is ~ 4', therefore all 
of the clouds are beam diluted and it is necessary to consider their beam-filling fractions 
to understand their relative intensities. All else being equal, larger clouds would appear 
brighter in Fig. ^. The peak of the CO ( J = 4 ^ 3) map is N159W, due to a combination of 
intense emission and beam filling. Although there is an extension of (J = 4 — i> 3) emission 
in the southward direction, a peak at the position of N159S is most noticeably absent. 
Because N159W and N159S have similar sizes and intensities in the CO (J = 1 ^ 0) 
transition, this implies that the CO (J = 4 — >^ 3)/(J = 1 — *• 0) ratio is much smaller for 
N159S. 

Conversely, there is a hint of an emission peak at the position of N160, a small 
cloud in ^^CO (J = 1 — > 0). The presence of this peak suggests that the CO 
( J = 4 ^ 3)/( J = 1 ^ 0) ratio for N160 is large. These results for the ratios are confirmed 
by convolving the CO (J = 1 ^ 0) map to the resolution of the (J = 4 ^ 3) data 
(Table H). The CO (J = 4 ^ 3)/(J = 1^0) ratio strongly depends on temperature and 
density, but the analysis performed by Heikkila et al. (1999) shows that both clouds have 
volume densities above the (J = 4 ^ 3) transition critical density. Therefore the deficit of 
( J = 4 ^ 3) emission in N159S is mostly due to the much lower temperature of the southern 
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cloud. In the optically thick, LTE limit the observed CO ( J = 4 ^ 3)/( J = 1 — 0) ratio 
corresponds to T^in ~ 8 K. This low temperature is probably due to the absence of star 
formation activity in N159S. 

Figure ^ reveals the sources of heating and photodissociating radiation in the 
N159/N160 complex. The optical N159 nebula is centered between the N159E and N159W 
clouds. Portions of these clouds are in the foreground of the nebula and show themselves 
as obscurations against the bright background. The N160 cloud, however, appears to be 
behind its nebula. Part of the N160 nebulosity fills in the gap between the N160 and 
N160E molecular clouds, and is neatly delineated by them. N159S is seen here only as a 
subtle obscuration against the field stars, with no associated nebulosity. While all the other 
clouds have apparent embedded sources visible as peaks in the 60 /xm IRAS HIRES picture, 
N159S appears completely devoid of such sources. This lack of 60 fim emission confirms the 
complete absence of any massive star-forming activity in the southern cloud. 

3.1.2. The Nature of the Extended CO Emission 

To investigate the faint extended emission apparent in the ^^CO (J = 2 ^ 1) map, 
two small subregions at the edges of the complex (squares in Fig. |l|) were mapped 
simultaneously in ^^CO ( J = 2 ^ 1) and (J = 1 ^ 0). This permitted us to obtain reliable 
line ratios for these two transitions. The average spectra for these two regions are shown in 
Fig. |. The ^^CO (J = 2 ^ 1)/(J = 1 ^ 0) ratio is ~ 3 for both positions; about a factor 
of 2 — 4 larger than for the molecular peaks (cf.. Table 0). For thermalized, optically thick 
CO, we expect a line ratio close to one. Larger (J = 2 — ^ 1)/(J = 1 0) ratios can be 
produced by essentially only four mechanisms: 1) self or foreground-absorbed (J = 1 — 0) 
emission, 2) optically thick CO gas with temperature gradients, 3) different beam filling 
fraction for the (J = 1 — > 0) and (J = 2 ^ 1) transitions, or 4) thermalized but optically 
thin CO emission. 

In the first of these cases the lower transition is absorbed by either a cold cloud along 
the line of sight at similar velocities, or colder outer layers of the same cloud. If most 
of the CO in the absorbing gas is in the ground state (i.e., the intervening cloud is very 
cold) the (J = 2 ^ 1) transition will not suffer this effect, thereby artificially raising the 
( J = 2 ^ !)/( J = 1 ^ 0) integrated intensity ratio. In case 2) the r = 1 surface arises 
at different places along the line of sight for the two transitions. The observed line ratio 
will be the ratio of temperatures of the regions where the r = 1 surface occurs. If there are 
temperature gradients in the cloud, then this ratio can in principle take any value. In case 
3), because the opacity grows faster with Nco in the ( J = 2 ^ 1) than in the (J = 1 — 0) 
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transition (Eq. it is possible for the clumps to appear larger in the higher transition 
and consequently fill more of the beam. For this effect to be considerable it requires small 
and warm CO clumps (cf., Appendix ^). In case 4) the ratio of both transitions will be 
in the range — 4, depending on the temperature of the emitting CO gas (cf., Fig. |^). 
It can be shown that, in general, the ratio of integrated brightness temperatures for two 
consecutive rotational transitions in the optically thin, thermalized limit is 



'.(J+l)^J /'J+l\^ _ h.,o{J+l) 



Tj^{j-i) \ J 



e 



where J is the rotational quantum number, Tex is the excitation temperature, and 
huio/k ~ 5.5 K for ^^CO. For the transitions considered, a line ratio of ~ 3 implies 
excitation temperatures T^x ^ 40 K. 

Which one of these possibilities is occurring in regions Rl and R2? The (J = 1 ^ 0) 
spectrum for R2 (Fig. ^ shows no clear evidence for self-absorption. Also, self-absorption 
tends to be a rather localized phenomenon: it is difficult to imagine it happening for two 
completely unrelated regions like Rl and R2. Thus we think it is an unlikely explanation. 
To decide among the remaining possibilities it would be ideal to be able to estimate the CO 
column density. The ^^CO observations, however, are not sufficiently sensitive to provide 
useful limits for the ^^CO/^'^CO ratio and subsequently determine if the emission is optically 
thick. 

Because some of the possibilities we are discussing require warm gas, it is important 
to consider the sources of heating in these two regions. Region Rl is located between the 
N160 and N159 nebulae, and consequently its temperature can easily be > 40 K. Region R2 
is found, however, in a quiescent region about 4' (~ 60 pc projected distance) away from 
the center of the N159 nebula. To raise the gas temperature of region R2 to the ~ 40 K 
needed to explain our ratios with optically thin CO emission, it is necessary to have heating 
sources. Surface temperature calculations of the PDR show that only a modest radiation 
field, Xuv, is required to produce that temperature (xuv~ 10, Kaufman et al. 1999). 
Eastward of N159S there are a few very faint and inconspicuous H II regions (DEM 272, 
277 and 279, see Fig. |l|; Davies et al. 1976), and traces of low level 60 /im emission (Fig. 

This suggests that photons, possibly from these H II regions or from their much brighter 
northern cousins, find their way there to heat the dust and elevate the temperature of the 



diffuse gas. Calculations in § p.2| show that the N159 nebula is probably bright enough to 
provide the Xuv~ 10 field needed to raise the temperature of the gas to 40 K. This radiation 
is not, however, intense enough to heat the core of N159S and make it a strong 60 fim or 
CO {J = 4 —> 3) emitter (Fig. It is certainly not enough to produce any measurable 158 
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/iin [C II] emission, a transition that requires T ~ 92 K to be excited (Fig. 

Given the fact that regions Rl and R2 are at the edges of the molecular cloud complex, 
we find the case for optically thin CO emission more compelling than the alternative 
explanations. It is, however, difficult to choose among possibilities 2), 3), and 4) with the 
available data. Nevertheless, we can make some specific predictions for the last two cases. 
In case 3) (larger beam-filling fraction for the (J = 2 ^ 1) transition), simple geometric 
arguments show that very small CO clumps are required (Appendix 0). These clumps are 
so small that they have only r ~ 2 in the (J = 1 ^ 0) transition and fill less than 1% of 
the beam. Future measurements of the (J = 3 ^ 2) and (J = 4 — >■ 3) transitions in these 
regions may help distinguish between case 3) and optically thin emission (case 4). Because 
the ratio of opacities between the (J = 3 ^ 2) and (J = 2 — ^ 1) transitions of CO is 
2.25 for infinite temperature (Eq. |ip, compared with 4 for r(J = 2 ^ l)/r(J =1^0), 
the increase in beam filling fraction going from the (J = 2 ^ 1) to the (J = 3 ^ 2) 
transition will be only modest. Thus, assuming uniform density clumps in LTE we expect 
peak intensities Tmb ~ 0.18 and 0.30 K in the ^^CO (J = 3 ^ 2) transition for regions Rl 
and R2 respectively, only 20% brighter than the observed (J = 2 1). For clumps with 
density increasing toward the center the difference between the beam filling fraction for 
both transitions, and consequently the increase in brightness temperature, will be even 
smaller. 

Concerning case 4), optically thin CO, the observed (J = 2 — * 1) intensities require 
CO column densities of Nco ~ 2 x 10^^ cm~^ assuming optically thin LTE emission at 
Te^. = 40 K. A prediction of this model is that ^^CO (J = 3 2) and (J = 4 ^ 3) 
emission should be readily observable with peak intensities in a ~ 1' beam of T^b ~ 0.25 
and 0.45 K for regions Rl and R2 respectively. Because of the low angular resolution of 
our ^^CO (J = 4 ^ 3) data, however, we are not able to cleanly separate these regions 
from the nearby peaks that dominate the emission and therefore cannot test these models. 
Nevertheless, the fact that these regions are relatively bright in [C I] ([ ^3.2.2| ), together with 
the results of column density calculations (^ |3.3D , strongly suggest that they are part of a 
translucent, mostly photodissociated envelope. 

3.2. Carbon in the Gas Phase of the ISM 

Here we analyze the AST/RO [C I] (^Pi ^^Pq) map in conjunction with the available 
[C II] and CO data (Israel et al. 1997; Johansson et al. 1998). Figure ^ reveals the 
distribution of the three dominant forms of carbon in the gas phase of the ISM. The striking 
features of this map are the complex interplay between CO, C° and C"*", and the fact that 
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[C I] peaks in the southern cloud, a quiescent region entirely devoid of strong UV sources. 

Throughout the map there is an overall anticorrelation between the distributions of 
[C II] and [C I]. Only in N159W do the three species peak at approximately the same place, 
while in all the other clouds only two of the species show intense emission. The northern 
regions, with abundant star formation activity, are bright in [C II] but very dim in [C I] . 
The opposite is true for N159S. This cloud, with no massive star formation activity, is the 
peak of [C I] in the whole complex. It is also the region with the largest [C I] /CO intensity 
ratio (cf.. Table The [C I] emission in N159E is relatively faint, and the CO, C° and C"*" 
there appear to occupy adjacent and partially overlapping regions. Finally, in N160 a [C I] 
hole is filled in by a bright lobe of [C II] , with a very faint ridge of [C I] emission overlapping 
with the CO. 

The distribution of the different forms of gas phase carbon in N159/N160 may be 
affected by peculiarities of the region. For example, the central [C II] peak (N159-M, Israel 
et al. 1996), which has no FIR counterpart, is very close to the position of LMC-Xl. It 
may be associated with that X-ray source, which is located at offsets ~ [— 1', +1'] (Fig- 
LMC-Xl, one of the strongest X-ray sources in the LMC, is a black-hole candidate with an 
07 III companion star (Cowley et al. 1995; Schmidtke, Ponder, & Cowley 1999). It has 
been suggested that X-ray radiation has considerable effect on the chemistry of the ISM 
(e.g., Lepp & Dalgarno 1996), and it is possibly a way to produce C"*" and C°. Dissociation 
of CO in shocks may be another way of generating these species. Approximately 2' east of 
LMC-Xl, Chu et al. (1997) have identified a supernova remnant, SNR 0540-697, based on 
optical spectra and X-ray data. This remnant is expanding at ~ ±150 km s^^ and overlaps 
with most of the nebulosity NW of N159E (Fig. |^). We see no evidence in our data for [C I] 
emission at these velocities, and the ratio maps (Fig. |^) show no peculiar enhancement 
of the I[ci]/Ico and I[cii]/Ico ratios at the position of the remnant. Thus, shock-induced 
dissociation of CO does not appear to be an important mechanism producing [C II] or [C I] 
in this area. 

On the spatial scale of these observations (1' ~ 15 pc) we do not expect to be resolving 
the PDRs into their three separate C"*", C° and CO regions. This has been accomplished 
only in a few Galactic sources observable with much greater spatial resolution (e.g., the 
Orion bar). Accordingly, we expect to observe coextensive [C II], [C I] and CO emission for 
the molecular peaks. This is not necessarily the case for the translucent medium. In diffuse 
gas and in the presence of strong UV sources most carbon will be ionized. Thus C"*" will 
dominate the emission, forming [C II] regions akin to Stromgren spheres around the ionizing 
sources. This is apparently happening in the northern portion of the map where bright 
lobes of [C II] emission east and west of N160 fill in holes in the [C I] distribution. These 
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[C II] peaks and [C I] holes are unequivocally associated with H II regions (Fig. Most 
of their [C II] emission may arise from C"*" mixed with H II inside the Stromgren sphere, 
coUisionally excited by electrons. 

In principle, over a small range of extinction (Av~ 0.3 — 1) we could have clouds where 
C+ and C'' are dominant, with little or no CO emission. That may be occurring between 
N160 and N159W, where bright [C II] and [C I] overlap in a region that shows little CO 
emission. One must be careful, however, when interpreting and combining data with very 
different (and, for the [C I], low) angular resolutions. Figure ^ shows the three transitions 
convolved to a common resolution (~ 4.3'). In this map there are only two peaks in the 
[C II] distribution (N160 and N159W), and the [C II], [C I] and CO maxima near N159W 
are displaced by about 1' from each other in a way that resembles a PDR C^/C^/CO 
transition. This structure, however, is not a PDR: if these clouds were moved to Orion (0.5 
kpc away and perhaps the best example of a resolved PDR) the distance between the [C II] 
and the CO peaks would span 2 degrees in the sky. What we are observing are large scale 
excitation and chemical gradients. The extended [C II] emission is heavily weighted in the 
convolved data and pulls the overall maximum northward, also diluting the peak coincident 
with N159W. In much the same manner, the (J = 1 — > 0) peak is pulled southward by the 
extended CO emission. Notice that the CO maximum near N159S moves westward for the 
same reason. 



3.2.1. The Neutral Carbon Emission from N159S 

The fact that the peak of [C I] emission for the entire complex is a quiescent region is 
unexpected if C° has, chiefly, a PDR origin and thus requires UV photons to be produced. 
To review the evidence: N159S is a dark cloud with weak CO ( J = 4 ^ 3) emission (hence 
at a low temperature), with no conspicuous heating sources apparent in the optical or the 
FIR (hence little or no star formation), and with very faint [C II] emission (hence no UV 
sources). Nevertheless, according to the analysis by Heikkila et al. (1999) it is the cloud 
with the largest column density in the complex {Nh2 ^ 1-7 x 10^^ cm~^, Av~ 4.5 compared 
to Nh^ ~ 1.1 X 10^2 cm-2, Av~ 2.9 for N159W and N160), and it is also the brightest [C I] 
emitter, possessing the largest I[ci]/Ico ratio. This is certainly indirect evidence, but it 
suggests that most of the C° in this cloud originates not by UV photodissociation of CO in 
the PDR, but inside the CO cores. 

Recent modeling results indicate, however, that the I[ci]/Ico ratio is very insensitive 
to the radiation field (Kaufman et al. 1999). Therefore only a small amount of UV 
radiation is necessary to explain the observed ratio in the context of a PDR. This can 
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be understood in the following way: in Ay space the extent of the [C I] emitting region 
in a homogeneous one- dimensional calculation, and consequently the C° column density, 
is relatively insensitive to the input radiation field. The depth at which the C"*'/C°/CO 
transition occurs, however, increases for larger Xuv Because ^^CO becomes optically thick 
soon after this transition, neither the [C I] nor the CO emerging intensities are a strong 
function of the radiation field at the surface of the cloud. All the molecular peaks in this 
region feature optically thick ^^CO, based on the ^^CO/^^CO intensity ratios listed in Table 
I (^2c/i3c~ 50 for the LMC, Johansson et al. 1994). 

The previous reasoning applies to homogeneous clouds with plane-parallel (i.e., 
one-dimensional) geometry. If molecular clouds are in fact clumpy, they can be modeled as 
an ensemble of spherical cloudlets. In this scenario the UV field should have a dramatic 
effect on the C^, C° and CO column densities and intensity ratios. In a clump subjected to 
an increasing UV field, the CO emitting region can be pushed towards the center only as 
far as the radius of the clump, after which the entire clump photodissociates. This scenario 
has been modeled in detail by Bolatto et al. (1999). Following those calculations we expect 
an inverse dependence between Xm and CO column density, and a modest increase in the 
[C I]/CO ratio for increasing Xm (growing by ~ 25% per order of magnitude in Xuv)- This 
model, however, does not include an interclump medium. Modeling of clumpy PDRs as a 
collection of high-extinction (Av~ 100), dense regions immersed in a diffuse interclump 
gas exposed to high Xm (Meixner & Tielens 1993, 1995) suggests that most of the CO 
(J = 1 — s> 0) and [C I] emission comes from the low density interclump medium. Thus, at 
low UV fields and low clump filling fractions, and when considering interclump dominated 
lines, a clumpy PDR may behave very similarly to a homogeneous plane-parallel model 
of interclump density. Similarly, at high clump filling fractions the contribution from the 
interclump medium would be unimportant. 

We will now try to estimate the radiation field in the N159S region. This cloud is 
most likely illuminated by the bright northern H II regions, although there may be a 
small contribution from the much smaller nearby regions mentioned in jj 3.1.2| (Figs. 0, 



Estimates for the luminosities and possible exciting sources of the FIR peaks in the 
N159/N160 complex, based on IRAS HIRES data, are given in Table |[ These estimates are 
in reasonably good agreement with the KAO results by Israel et al. (1996). An equivalent 
spectral type and multiplicity is assigned to each flux, by assuming that all of the starlight 
is trapped by the dust and reradiated in the FIR. 

In order to use these results to estimate the UV radiation field far away from these 
sources we will assume that the radiation intercepted by dust is only a fraction / of the 
total luminosity, and therefore a fraction (1 — /) escapes the star-forming region and the 
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enshrouding dust. Furthermore, we will assume that about 50% of the luminosity of the 
exciting stars is emitted in the range 912 — 2066A (typical for late 0/early B-type stars). 
Therefore the UV radiation field at a distance d is 



(1 - f)LFIR 

where xo is the standard normalization factor (approximately the UV radiation field in the 
vicinity of the Sun, xo = 1-6 x 10~^ erg s~^ cm~^; Habing 1968). Using this formula we can 
compute Xuv for N159S, assuming that it is illuminated by the nearby N159 H II region. At 
the projected distance c? ~ 90 pc, and assuming that 90% of the UV escapes the surrounding 
dust, we obtain Xm — 60. If only 50% of the UV escapes then the radiation field would be 
Xuv — 6. By comparison, radiation field estimates by Israel et al. (1996) for the northern 
clouds N159E-W and N160 near H II regions range from 300 — 600 (cf.. Table p. Our Xm 
estimate for N159S is probably too large, since: 1) we use the projected distance, and 2) we 
assume no dust extinction between the northern H II regions and N159S. 

Another way of computing the luminosity of the UV sources in this region is to use 
the radio continuum information. The 21 cm radio continuum data (Figure ^ poses 
constraints similar to the FIR on the luminosity of the stars that excite the H II regions. 
In Table ^ we have calculated the spectral type and multiplicity, assuming optically thin 
free-free emission and following Jackson & Kraemer (1999) analysis. These estimates are 
generally lower than those based on the FIR, mostly because spectral types later than 06 
are still very luminous (adding to the FIR luminosity) but contribute little to the radio 
continuum. Recent ISOCAM observations (Comeron & Claes 1998), for example, reveal 
three strong 15 /xm peaks in N159 which feature very faint radio continuum, and are thus 
attributed to ultra-compact H II regions which are optically thick at 21 cm. These are 
examples of sources that would add to the FIR luminosity but be invisible at 1.42 GHz. 
Nevertheless, both the radio continuum and the FIR measurements agree within a factor 
of ~ 2 {Lpjji ~ 2Ljic), strongly suggesting that most UV is actually intercepted by the 
surrounding dust. Taking into account the various estimates and caveats discussed in the 
previous paragraph, we conclude that Xuv(N159S);^ 10. Such a low UV field is consistent 
with the faint [C II] emission from N159S, as well as its low I[cii]/Ico ratio. 

Even orders of magnitude in Xmy however, have little impact on the value of the 
I[ci]/Ico ratio in standard PDR plane-parallel models. According to the calculations by 
Kaufman et al. (1999), the [C I]/CO intensity ratio observed in N159S can be attributed to 
gas with n ~ 10^ cm~^ at Xuv~l — 10^. Note that metallicity has only very small effects on 
plane-parallel model calculations, for reasons essentially similar to those discussed for Xm- 
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Thus, although the aforementioned calculations were performed for Galactic sources their 
results can be safely applied to the LMC The coincidence between this density and that 
obtained by Heikkila et al. (1999) for N159S appears satisfying, until we take into account 
the different scale of the measurements: the 4' (60 pc) region over which I[ci]/Ico— 0.15 is 
^ 25 times larger than the SEST ~ 48" beam used to carry out the multiline excitation 
analysis. This is an extremely high density for such a large region. Invoking beam filling 
fraction arguments avoids this problem, by assuming that only small and dense regions 
within the beam are dominating the CO ( J = 1 — > 0) and [C I] emission. Comparison of 
the measured [C I] intensity with the model results suggests ~ 10% beam filling fraction. 
Nevertheless, if such a high density is characteristic of the gas producing the [C I] emission, 
it rules out the interclump medium as the main reservoir of neutral carbon. Meixner & 
Tielens' (1995) modeling of two-phase (clump + interclump) PDRs shows that I[ci]/Ico~ 10 
can be obtained for a mixture of 5 x 10^ cm~^ Av~ 100 clumps immersed in a 3 x 10'^ cm~'^ 
diffuse component, with 20% volume filling fraction for the clumps. These computations, 
however, were carried out for Xuv~ 5 x 10^ and no similar calculations are available for UV 
fields closer to that of N159S. 

At the beginning of this section we pointed out a series of intriguing facts about N159S 
that may suggest a non-PDR origin for most of its neutral carbon. Despite them, however, 
we conclude that the observed [C II], [C I], and CO intensities in N159S are consistent with 
standard PDR theory assuming n ~ 10^ cm~^ and Xuv~ 10. These values of the density 
and UV field agree with previous excitation analysis and with the best estimates of Xuv in 
N159S, assuming it is illuminated by the northern complex of H II regions. 



3.2.2. The I[ci]/Ico Ratio in the Complex 

How are the relative intensities of the different forms of carbon affected by the local 
conditions? In the previous sections we have seen that [C I] is dim in the northern portion 
of the map, where there is active emassive star formation, and bright in N159S. In this 
section we will compare the I[ci]/Ico and I[cii]/Ico ratios throughout this region, and look 
for the effects of radiation fields and metallicity. 

Figure |^ shows the I[ci]/Ico and I[cii]/Ico ratios mapped for the entire complex, 
with the data convolved to a common 4.3' resolution. Table ^ shows our I[ci]/Ico results 
tabulated for the regions of interest. Notice that Table ^ uses a 4' beam size. To convert 
from intensity ratios given in K km s"^ to ratios in erg s~^ cm~^ sr~^ multiply by ~ 78 for 
I[ci]/Ico- Since N160 is such a weak [C I] source, its ratio is dominated by the surrounding 
extended [C I] emission and therefore the ratio increases noticeably when measured in a 
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larger beam. For N159S the angular resolution of the measurement is unimportant, and 
the ratio remains mostly unchanged. Therefore the high I[ci]/Ico ratio at N159S is not an 
artifact caused by the nearby emission. The first measurement of the [C I] /CO intensity 
ratio in N159 was carried out by Stark et al. (1997a), who obtained I[ci]/Ico— 0.26 
(intensities in K km s~^). This determination was based on a single spectrum taken towards 
the nominal position of N159W, and assumed an unresolved structure for the CO from this 
cloud. We find that, albeit in a larger beam, I[ci]/Ico— 0.12 for N159W. This difference is 
largely due to an underestimate of the CO intensity in the aforementioned paper. 

Two regions of the map feature high [C I] /CO intensity ratios: the northern region, 
where there I[cii]/Ico also peaks, and the southeast corner, where there is extended low-level 
CO (J = 2 ^ 1) but little (J = 1 — > 0) emission and no detectable [C II] emission. We 
think that in both cases we are seeing translucent, mostly photodissociated gas, with the 
important difference that the abundant UV radiation in the northern region raises the 
temperature of the PDRs over the 91.3 K required to excite the 158 fim [C II] transition 
(see § p73|) . There is the possibility that the southeastern extension of the N159S cloud in 
[C I] is an artifact of the sampling, since the map is missing a few spectra there. Because 
of the heavy oversampling of the beam, however, we think that this extension is probably 
real. It appears also that there is a modest increase in the I[ci]/Ico ratio at the edges of the 
complex, and perhaps in the immediate vicinity of LMC-Xl. The lowest ratios are found 
south of N159W. 

How does the I[ci]/Ico ratio in this complex compare with the I[cii]/Ico ratio? Since 
both C^ and C° are produced by the action of UV photons in the PDR, the naive 
expectation is that [C II] and [C I] should have a similar distribution in unresolved PDRs. 
Excitation differences, important for very low radiation fields (Xuv< 10), will become 
negligible for Xuv> 100 as the PDR reaches temperatures well over 100 K (Kaufman et al. 
1999). As discussed at the beginning of §3.2| there are other reasons why we may not expect 
[C II] and [C I] to be coextensive in the diffuse translucent gas, namely the formation of 
C^ "Stromgren spheres" near UV sources. In these regions carbon ionization and not dust 
absorption is the dominant process that removes UV photons, and consequently all the 
carbon is in the form of C'^. In the molecular material, however, UV photons should be 
predominantly removed by dust grains, leaving a fraction of the carbon in the form of C^. 
Consequently, if the structure of the ISM is clumpy and therefore allows the UV photons to 
penetrate deep into the clouds, the expectation is for [C I] and [C II] to be coextensive and 
in many ways behave similarly. In particular, modeling of clumpy PDRs suggest that both 
the I[cii]/Ico and I[ci]/Ico ratio should be enhanced in low metallicity systems (Bolatto et 
al. 1999). This is mostly due to the rapidly dwindling sizes of the CO cores of the clumps 
with decreasing metallicity and dust-to-gas ratio. 
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The [C II] /CO intensity ratio has been studied in a variety of systems (e.g., Stacey 
et al. 1991), and the highest ratios are associated with low metalhcity environments. As 
we pointed out in §0, Mochizuki et al. (1994) found a considerable enhancement of the 
global I[cii]/Ico ratio of the LMC (I[cii]/Ico~ 23, 000) over that of the Milky Way. Average 
ratios for Galactic objects are I[cii]/Ico~ 1,300 for GMCs, and ~ 4,400 for H II regions 
(Stacey et al. 1991). Observations by Israel et al. (1996) of N159/N160 (cf.. Table ^ 
show that the three northern molecular peaks (N159E, N159W and N160) have I[cii]/Ico 
ratios 1 — 5 times larger than comparable Galactic regions while their Ifir/Ico ratios are 
3 — 40 times lower, suggesting a much larger abundance of Measurements of the 30 
Doradus region (Poglitsch et al. 1995) yielded a ratio I[cii]/Ico~ 65,000 for its molecular 
peaks. Madden et al. (1997) studied the global [C II] emission from the low metallicity 
dwarf galaxy IC 10 which has Ziao ~ Zmw/4 (where we have used the oxygen abundance 
in Orion as representative of the Galaxy, 12 + log[0/H] ^ 8.75; Lequeux et al. 1979). They 
found ratios in the range I[cii]/Ico~14,000 — 87,000 for various regions. 

Unlike the widely varying I[cii]/Ico ratio, the I[ci]/Ico ratio for the molecular peaks of 
N159/N160 is surprisingly uniform. The values range only between I[ci]/Ico~ 7 to ~ 12. 
These are very similar to the typical ratio in the Milky Way where I[ci]/Ico~ 13 (intensities 
in erg s~^ cm~^ sr~^). This is not an isolated result. Wilson (1997) observed the [C I] 
emission from several clouds at different galactocentric distances in M 33, a spiral galaxy 
with a well-known metallicity gradient. She obtained ratios of I[ci]/Ico~ 3 — 14 and found 
no obvious trend with inferred metallicity. Recently, Bolatto et al. (2000) studied the 
I[ci]/Ico ratio in the molecular cloud complex associated with the brightest star- forming 
region of IC 10. They obtained an average ratio I[ci]/Ico~ 18, only slightly larger than 
the average Galactic ratio. Finally, Gerin & Phillips' (2000) recent study of atomic carbon 
in a variety of galaxies finds some dispersion in the I[ci]/Ico intensity ratio, but no clear 
segregation with galaxy type. Their sample of galaxies has an average ratio I[ci]/Ico~ 16, 
with most ratios found the in the interval I[ci]/Ico ~8 — 32. 

Overall, if there is a trend for I[ci]/Ico with metallicity it is certainly much weaker than 
the one observed for I[cii]/Ico- A plausible explanation for the constancy of the [C I] /CO 
intensity ratio is a non-FDR origin for a fraction of the C°. While it seems unequivocal 
that in the diffuse gas most neutral carbon is associated with FDRs, it is possible that a 
different mechanism dominates its production inside molecular cloud cores. If C° is indeed 
produced in these cores by processes unrelated to photochemistry, it would simply explain 
the close observed association between the [C I] and CO line intensities. 
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3.3. The Column Density in the Extended Envelope 

The entire molecular cloud complex appears to be surrounded by an extended envelope 
visible in CO (J = 2 — > 1), and possibly [C I] and [C II]. For example, the extended and 
relatively bright [C I] emission east of N159S and south of N159E correlates very well with 
the faint, extended CO ( J = 2 ^ 1) in the same region, as does the tongue of [C I] spreading 



north of the N159 nebula. Assuming, as was discussed in § p . 1 . 2| , that the CO emission in 
this region is optically thin, then we appear to be seeing translucent gas emitting in [C I]. 
Translucent clouds are clouds with visual extinction Ay^ 1, where an important fraction 
of the carbon is in the form of C"*" and C° (e.g., Ingalls et al. 1997). While the regions 
surrounding Rl are strong [C II] emitters, region R2 and its environs are very faint in [C II]. 
Undetected by Israel et al. (1996), the upper limit for the integrated intensity from R2 is 
I[cii]~ 6.8 X 10~^ erg s~^ cm~^ sr~^. The limit for the corresponding C"*" column density, 
derived for the high-temperature, high-density case (that is, the minimum column density 
corresponding to the intensity limit) is Nc+ ^ 4 x 10^^ cm~^, about 20 times larger than 



the corresponding column density of CO derived in §3.1.2. Singly ionized carbon could well 



be the dominant form of carbon in this region yet remain undetected. 

Assuming that most of the carbon is in the form of C+, we can estimate the hydrogen 
column density along these lines of sight. For Rl, l[cii] ~ 1.4 x 10^^ erg s~^ cm~^ sr^^, 
therefore Nc+ ~ 8 x 10^^ cm~^. Assuming an LMC gas phase carbon abundance per H 
nucleus xc ~ 7 x 10"^ (Dufour 1984), this yields Nh ~ 1.2 x 10^^ cm"^, or Ay ~ 0.15 (using 
the conversion factor lAy ~ 8 x 10^^ cm~^ which is four times that of the Galaxy). In 
the case of R2, where no [C II] is detected. Fig. ^ shows that there is some very low level 
emission that becomes statistically significant only after spatial smoothing. We will assume 
I[cii]~ 3 X 10~^ erg s~^ cm~^ sr~^ (about half the sensitivity limit). In the high-temperature, 
high-density limit this would result in Nh ~ 4 x 10^° cm~^, or Ay ~ 0.05. If we allow for a 
relatively low PDR temperature, T^in ~ 40 K, consistent with a low Xnv and the optically 
thin CO analysis, we find Ay ~ 0.15. It is important to point out that these are line-of-sight 
extinctions, and not necessarily equal to the Ay that extincts the UV field and enters the 
PDR calculation. For a completely edge-on PDR, for example, these two extinctions would 
be measured along orthogonal directions and thus be unrelated. If this gas fills the beam 
and has a density near the critical density of [C II] {ricr ~ 3 x 10^ cm~^), as we implicitly 
assumed in the calculations, then we are seeing a ~ 0.1 pc thick layer of gas. We find this 
sheet-like geometry, with a sheet thickness of only 1/100 of the extent in the plane of the 
sky, uncomfortable although not impossible. 

How do these results compare with column density predictions based on the H I data? 
We assume a spin temperature of Ts„ ~ 110 — 180 K based on the H I absorption observed 
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towards the N159 continuum sources (Spitzer 1978). We can then compute atomic hydrogen 
cohimn densities of Njj ~ 4.2 x 10^^ cm~^ and Nh ~ 4.5 x 10^^ cm~^ for regions Rl and 
R2 respectively. These column densities are ~ 3 times larger than those derived from I[cii]- 
The gas is apparently hot enough to populate the upper fine structure level of C"*" and 
excite the 158 fira [C II] transition (recall hv/k = 91.3 K for [C II]), thus this discrepancy is 
probably not due to temperature effects. 

This method, however, is not very precise since it heavily relies on similar excitation 
conditions for the H I seen in absorption and emission, and this assumption is compromised 
by the large scales involved in the map (recall l'~ 15 pc). Other possible causes for the 
apparent discrepancy in Nh are: 1) a deficiency of C+ in the atomic gas (only 30% of 
the C is photoionized) , 2) the volume density of the atomic gas is low, and consequently 
the excitation of [C II] is subthermal {n < ncr{[C ll])), or 3) there are density fluctuations 
within the beam, and the [C II] emission is dominated by the fraction of the C"^ column 
that is thermahzed {n > ncr{[C II])). 

The flrst possibility appears to be very improbable, since unless the material is well 
shielded from the UV most of the carbon should be in the form of C"*" rather than C°. 
Unfortunately the large beam in our [C I] observations makes impossible to separate 
cleanly these regions from the molecular peaks, thus precluding us from obtaining a reliable 
estimate for the column density of in Rl and R2. Subthermal excitation of requires 
the density of H I to be below the critical density of the 158 fim. transition, ricr ~ 3000 
cm~^. The intensity l[cii] of C"*" coUisionally excited by H atoms can be computed using 
(e.g.. Madden et al. 1997): 



I[cii] = 2.35 X 10-21 Nc+ 



2ex-p{-hp/kT) 



1 + 2exp{—hiy/kT) + ricr/n 



(3) 



where hu/k = 91.3 K. Consequently, a discrepancy of a factor of 3 would require n ~ 500 
cm"^, well below the critical density of [C II]. Notice that this density would bring the 
physical thickness of the [C II] emitting layer to ~ 3 pc. The third possibility (density 
fluctuations within the beam, akin to clumping) would make the regions where n > ricr 
dominate the [C II] emission (about 1/3 of the total column), while in most of the gas 
n <^ Tier and the [C II] excitation is subthermal. This ought to be happening at some 
level, because CO ( J = 2 — > 1) is present throughout the region and its excitation requires 
n ~ 10^ cm"^. This density is probably too large to be the average volume density in 
the envelope, thus most of the CO and perhaps a large fraction of the [C II] emission is 
originating in clumps within the envelope. 



4. Summary and Conclusions 



We have discussed new ^^CO (J = 1 ^ 0), ^^CO (J = 2 -> 1), ^^CO (J = 4 ^ 3), 
and [C I] (^Pi ^^Pq) emission line maps of the N159/N160 molecular cloud complex of the 
Large Magellanic Cloud. 

The ^^CO(J = 2^1) map shows extended faint emission previously undetected in 
the {J = 1 —> 0) transition. Further analysis of the {J = 2 —>■ 1)/{J = 1 ^ 0) intensity 
ratio in two selected subregions (Rl and R2) shows that this ratio is ;^ 3 for the extendend 
low level emission. This high value indicates optically thin CO emission from warm gas 
(Tkin > 40 K). The I[ci]/Ico ratio for Rl and R2 appears modestly enhanced with respect 
to the molecular peaks (Fig. ^. The I[cii]/Ico ratio is large for the northern region (Rl), 
but [C II] was not detected by Israel et al. (1996) towards region R2. Because of the low 
UV field heating the ISM in the southern region, this may be a temperature effect. The 
gas in the southern portion of the complex, far away from the massive star formation and 
bright H II regions, is too cold to excite the 158 /im transition (T < 92 K). The faint CO 
envelope appears to be translucent (Ay < 1), and column density calculations assuming 
that C^ is the dominant form of carbon confirm this result. The column density derived 
from the neutral hydrogen 21 cm data appears to be ~ 3 times larger than obtained 
from [C II]. This is probably caused by density enhancements (clumping) within the beam. 

The I^CO (J = 4 3) map shows that the (J = 4 ^ 3)/(J = 1 ^ 0) intensity 
ratio is ~ 4 times larger in the northernmost cloud (N160) than in the southern cloud 
(N159S). This agrees very well with indicators of star formation activity (Ifir, I[cii]; radio 
continuum), and is probably due to the much lower temperature of the quiescent southern 
cloud. Estimates of the radiation field incident on N159S suggest Xuv^ 10, consistent with 
a low temperature. The ratio of 100/60 fim continuum also indicates a low temperature. 

The [C I] ('^Pi ^'^Pq) map shows the [C I] intensity and the I[ci]/Ico ratio for the 
molecular concentrations peaking in N159S. The observed values are consistent with the 
PDR computations available for homogeneous, plane-parallel models, at n ~ 10^ cm~^ and 
Xnv ~ 10. Because in those models I[ci] is not sensitive to Xuv, however, the presence of 
intense [C I] emission in N159S is not a good test for a possible non-PDR origin for C° 
inside GMCs. There is pervasive [C I] emission throughout the region between N159EW 
and N160. This region also emits strongly in [C II], and features some of the largest I[ci]/Ico 
and I[cii]/Ico ratios in the complex. The radio continuum sources embedded in this gas 
appear as holes in the [C I] distribution, while they are peaks in the [C II]. The [C II] in 
these peaks probably originates inside the H II regions, mostly excited by collisions with 
electrons. 
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The [C I] /CO intensity ratios measured in the molecular peaks of the complex range 
between I[ci]/Ico 7 — 12 (intensities in erg s~^ cm~^ sr"-*-). This ratio is similar to the 
average ratio in the Milky Way (l[ci]/Ico~ 13) and to that measured in the low metallicity 
dwarf IC 10 (l[ci]/Ico ~ 18, Bolatto et al. 2000). The [C I] /CO intensity ratio for the 
translucent gas is somewhat larger, l[ci]/Ico ~ 20 — 30. While intense [C II] (^^3/2 — ^^Pi/2) 
emission and large I[cn]/Ico ratios appear to be unequivocally associated with massive star 
formation, [C I] (^Pi — *>^Po) emission and the I[ci]/Ico ratio shows no such clear association 
and has a much more complex behavior. The I[ci]/Ico ratio is more uniform than the 
I[cn]/Ico ratio throughout the complex, where spans values of I[ci]/Ico ~5 — 25 whereas 
I[cii]/Ico ranges between ~500 — 25,000 (i.e., a factor of 5 versus a factor of 50). This is 
partially due to the more stringent excitation conditions required by [C II]. 

We believe that our understanding of the ISM, especially the actively star-forming 
ISM, can be enormously advanced by further studies of the Magellanic Clouds. These 
galaxies, with their proximity, their active star formation, their low metallicities, and 
their unobscured lines of sight, present unique opportunities for detailed multiwavelength 
studies. Among the different regions in the Clouds, the N159/N160 complex is located at 
a privileged place, between the violent starburst of 30 Doradus and the quiescent southern 
CO ridge. Because of its location, its distinct environments, and the wealth of phenomena 
taking place within its bounds, the N159/N160 complex is one of the most interesting 
places in the Magellanic Clouds. Unfortunately, these galaxies are out of the reach of the 
large radioastronomy facilities in the Northern Hemisphere. The advent of some of the 
observatories now being planned, such as the Atacama Large Millimeter Array or a large 
South Pole telescope, will change that. The study of the ISM in these objects will vastly 
benefit from large submillimeter telescopes, array receivers, and millimeter/submillimeter 
interferometers located south of the equator. Future studies of N159/N160 should include 
the detailed characterization of the translucent molecular envelope, the dust properties, and 
the heating and cooling balance of the ISM near and away from star formation sites. Other 
observational programs could address the creation of molecular clouds by ram pressure 
compression, and their disruption by star formation activity. 
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A. Beam filling ratio for two transitions 



For this calculation we will assume spherical clumps of radius R, immersed in an 
isotropic UV field. A transition turns optically thick (i.e., its opacity, r, is unity) at a 
distance x from the surface of the clump. Therefore the radius of the clump at given 
transition (i.e., the distance from the center of the clump to the transition's r = 1 surface) is 
r = R — X. The ratio of beam filling fractions $ in two transition will then be proportional 
to the ratio of the projected areas 



$2 fR-X2^^ 



$1 \R-xi 



(Al) 



It can be shown easily that for gas in LTE the ratio of optical depths in the ( J = 1 — 0) 
and ( J = 2 — > 1) transitions will be 
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= 2 







hu21 


— e 


kTex 






\ 







(A2) 



where h and k are Planck's and Boltzmann's constants respectively, and Tex is the excitation 
temperature. For uniform density clumps the r = 1 surface occurs first for the transition 
with faster growing opacity, that is, 0:2/0:1 = ti/t2. Equation KTI thus has the solution 



R = xi- 



11 

T2 



(A3) 



In order to reproduce the CO ( J = 2 ^ l)/( J = 1 ^ 0) intensity ratios observed in 
regions Rl and R2, $2/$! ~ 3. Using T2/T1 ~ 3 (i.e., T^^ ~ 40 K according to Eq. [A2D 
yields R ~ 1.9 xi. Thus r(l — > 0) ~ 2 at the center of the clump. The average over the 
whole projected spherical clump is ~ 30% larger, or r(l ^ 0) ~ 2.5. 
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Fig. 1. — Carbon monoxide main beam integrated intensity maps acquired at SEST. The 
^^CO (J = 2 ^ 1) contours (left panel) are 2, 5, 10, 15, . . . , 55 K km s~^ {<Jrms ~ 0.6 
K km s~^) and the data have been convolved with a Gaussian profile to 33" resolution. The 
^^CO ( J = 1 0) contours (right panel) are 0.6 to 4.4 at 0.6 K km s~^intervals {arms ~ 0.2 
K km s~^) and the map has been convolved to 60" resolution. The individual pointing 
locations are shown by small dots. The boxes show the position of the two regions with 
low level extended emission also observed in ^^CO (J = 1 — > 0). The map intensities have 
been integrated over the velocity range 225 to 245 km s^^. The black star denotes the 
position of the black-hole candidate and strong X-ray source LMC-Xl. The black circles 
mark the positions of DEM L272 and L279, the two H II regions closest to N159S. The ^^CO 
( J = 1 ^ 0) map has been annotated with the nomenclature used in Table ^ 



2 0-2 
R.A. Offset (arcmin) 

^^CO(1^)/''^CO(1^0) 



0.05 0.075 



0.1 



0.125 0.15 



2 0-2 
R.A. Offset (arcmin) 

^^CO (2^1)/^C0(1^) 



0.25 0.5 0.75 



1.25 1.5 



Fig. 2. — Comparison of ^^CO ( J = 2 ^ 1) and ^^CO {J = 1 0) emission with published 
^^CO {J = 1 0) data (Johansson et al. 1998). All maps have been Gaussian smoothed 
to a common resolution of 60". The contour levels are the same in both maps (2, 5, 10, 
15, . . . , 40 K km s"^). Note that the ^^CO (J = 1 ^ 0) map was acquired by frequency 
switching, and required an elaborate baseline subtraction scheme. Therefore the lowest 
contour is very unreliable. (Left) Main beam integrated intensity ^^CO (J = 1 ^ 0) contour 
data (Johansson et al. 1998) overlaid on the ^'^CO/^^CO (J = 1 ^ 0) ratio image. Blue 
denotes low ^'^CO/^^CO ratios and therefore optically thinner CO emission, whereas red 
shows high ratios and optically thicker emission. (Right) Main beam integrated intensity 
^^CO (J = 2 — > 1) contour data overlaid on the ^^CO (J = 2 — > 1)/(J = 1 — > 0) ratio image. 
Blue signifies low (J = 2 ^ 1)/(J = 1 — > 0) ratios, and therefore possibly cooler (less dense) 
gas; red zones are warmer (denser). 
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R.A, Offset (arcmin) R.A. Offset (arcmin) 

Fig. 3. — Heating of the molecular cloud complex. The optical (left panel, blue band, SERC) 
and 60 /xm FIR (right panel, IRAS HIRES, 40 iterations) emission are shown here with ^^CO 
( J = 2 ^ 1) contours overlaid (alternate black and white, same levels as Fig. |l]). This picture 
illustrates the detailed structure of the nebulae and the intimate and complex association 
between nebulosity, obscuration lanes, and molecular clumps. The southern cloud appears 
to have no UV or heating sources associated with it. 
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Fig. 4. — Color image of the ^^CO (J = 4 — 3) integrated intensity map with ^^CO 
(J = 1 — > 0) data overlaid. The (J = 4 3) data has been smoothed to a 4' spatial 
resolution. Marked by crosses are the positions of the four main molecular peaks in the 
region. Their corresponding ^^CO ( J = 4 ^ 3) (black) and [C I] (blue) spectra are shown on 
the side (the [C I] spectra have been scaled up by a factor of two) For purposes of comparison, 
shown in red are the ^^CO (J = 1 ^ 0) spectra (Johansson et al. 1998) convolved to the 
same angular resolution (4' HPBW). The relative enhancement and deficit of (J = 4 — > 3) 
in the N160 (labeled N) and N159S (labeled S) clouds respectively are apparent. 
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Fig. 5. — Observed spectra of the ^^CO(J =1^0) and (J = 2 ^ 1) transitions for the 
faint extended molecular emission. These two regions are identified here by their offsets and 
drawn on Fig. |l|. The (J = 2 -^■ 1) transition (fine black line) produces the strongest lines 
in both regions while the ( J = 1 ^ 0) (wide gray line) is weaker by a factor ^ 3, suggesting 
optically thin emission. 
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Fig. 6. — Ratio of (J = 2 ^ 1) to (J = 1 ^ 0) integrated brightness temperature for 
CO emission in local thermodynamic equilibrium (LTE), as a function of the excitation 
temperature Tex- The black line corresponds to optically thick emission. The grey line is 
the plot of Eq. |T| for n = 1, valid under optically thin conditions. 



-34- 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



Fig. 7. — Distribution of the three dominant forms of gas phase carbon in N159/N160. The 
pseudo color image shows the AST/RO [C I] observations at an angular resolution of ~ 4.3' 
(sensitivity Icr ~ 0.1 K km s~^). (Left) In black contours is the CO (J = 2 ^ 1) map, 
at a resolution ~ 33" (contours 2,5,15,25,. . . ,55 K km s~^). In white contours is the [C II] 
map by Israel et al. (1996), at a resolution ~ 55" (contours start at 6.8 x 10~^ in steps 
of 6.8 X 10~^ erg s~^ cm~^ sr~^). (Right) The three transitions have been smoothed to a 
common resolution of 4.3'. We use here the CO (J = 1 ^ 0) map by Johansson et al. 
(1994) to avoid edge effects when convolving the CO. The CO intensity contours start at 2 
K km s~^ in steps of 2 K km s^^. The I[cii] contours start at 1 x 10~^ in steps of 2 x 10^^ 
erg s~^ cm~^ sr^-*^. 
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Fig. 8. — Maps of the I[ci]/Ico and I[cii]/Ico ratios throughout the N159/N160 complex. 
(Left) Image of the I[ci]/Ico integrated intensity ratio, with Ipi] and Ico (J = 1 — >• 0) in 
erg s~^ cm~^ sr~^. The CO (J = 2 — > 1) contours are overlaid to serve as a reference, but 
the (J = 1 — > 0) map from Johansson et al. (1994) was used to produce the ratio, thereby 
eliminating edge problems in the CO dataset. (Right) Image of the I[cii]/Ico integrated 
intensity ratio. Typical Galactic ratios are I[ci]/Ico~ 10 and I[cii]/Ico~ 4400. All maps 
have been convolved to the same ~ 4.3' angular resolution. 
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Fig. 9. — Maps of neutral carbon emission and 21 cm radio continuum. The radio continuum 

contours (white) start at 30 mJy and are logarithmically spaced by a factor of 1.4. The 
highest contour is at 620 mJy. The I[ci] contours (grayscale) start at K km s~^ and go to 
1.6 in steps of 0.1 K km s~^. 
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Table 1. Positions and Properties of the Regions of Interest in N159/N160 



Identifier 


R.A. 


Dec. 


R.A. offset 


Dec. offset 


n(H2) 


Nh2 


Xuv 


I[cii]/Ico ^ 




(B 1950) 


(B 1950) 


(arcmin) 


(arcmin) 


(cm-3) 


(cm- 2) 




N159W 


5''40™02!2 


-69°47'06" 


-1.4 


-0.1 


3 X 10^ 


1.1 X 10^2 


600—250 


5600 


N159E 


5>'40'"35?8 


-69°46'12" 


+1.5 


+0.8 








13000 


N159S 


40™ 2717 


-69°52'00" 


+0.8 


-5.0 


1 X 10-^ 


1.7 X 10^2 


60—6 


400 


N160 


40™ 06? 9 


-69°40'30" 


-1.0 


+6.5 


5 X 10^ 


1.1 X 10^2 


450—200 


34000 


N160E 


5h40™46?l 


_69°41'48" 


+2.8 


+5.2 










Rl 


5'^40™13?8 


-69°43'48" 


-0.4 


+3.2 




6 X 10^0 


200—50 




R2 


5h4om46!2 


-69°49'48" 


+2.4 


-2.8 




2— 6x10^0 


100—10 





''Prom the mm multiline excitation analysis by Heikkila et al. (1999). 
''Prom Israel et al. (1996). 
''Prom this work. 

''To convert to visual extinction, Ay, divide by 4 x 10^^. 
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Table 2. Measured Millimeter Line Parameters in N159/N160 



Identifier 


T-i a,b 


Vlsr 


FWHM 


1= / dv ^ 


I/Iio "■'^ 






(K) 


(km s-l) 


(km s-i) 


(K km s-i) 




(J = 1 


0) " 












N159W 




0.64±0.01 


238.1±0.1 


6.9±0.1 


4.67±0.04 


0.12±0.02 


N159E 




0.60±0.01 


234.6±0.1 


e.iito.i 


2.13±0.04 


0.07±0.01 


N159S 




0.57±0.01 


235.6±0.1 


5.7±0.1 


3.37±0.07 


O.lOitO.Ol 


N160 




0.56±0.01 


237.7±0.1 


4.6±0.2 


1.39±0.06 


0.12±0.02 


Rl 










O.OliO.lO 




R2 










0.09±0.11 




i^CO (J = 2 


^1)« 












N159W 




4.80±0.01 


237.7±0.1 


7.8±0.1 


39.12±0.06 


1.01±0.14 


N159E 




3.60±0.01 


234.6±0.1 


7.3±0.1 


27.42±0.06 


0.89±0.13 


N159S 




3.22±0.01 


234.9±0.1 


S.litO.l 


27.00±0.08 


0.78±0.11 


N160 




3.46±0.01 


237.7±0.1 


5.0±0.1 


18.24±0.10 


1.59±0.23 


Rl 




0.16±0.01 


238.0±0.7 


10.8±1.5 


1.75±0.20 


3.57±1.40 


R2 




0.25±0.01 


235.4±0.2 


7.6±0.4 


1.94±0.10 


2.81±0.64 


''Errors are 


statistical and iblcr. Systematic calibration uncertainty is la ' 


- ±10%. 



'^Values and errors derived from Gaussian fit. 

•^Ratio of transition intensity to that of ^^CO( J = 1 — > 0), convolved to the same beam, 
in K km s~^. 

''Errors are ztlu, and include 10% la calibration uncertainty in both lines. 
■^Measured in a HPBW=l' beam. 
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Table 3. Measured Submillimeter Line Parameters in N159/N160 



Identifier 




rp a,b 

(K) 


(km s-i) 


FWHM ^ 
(km s-i) 


1=/ T^tdv'' 
(K km s-i) 


I/Iio "•'^ 


I[cii]/Ico " 


12CO (J = 


















N159W 






0.95±0.01 


237.7±0.1 


8.6±0.1 


8.43±0.08 


0.58±0.09 




N159E 






0.42±0.01 


235.6±0.1 


10.2±0.2 


4.42±0.09 


0.48±0.07 




N159S 






0.40±0.01 


235.5±0.1 


7.7±0.2 


3.41±0.08 


0.27±0.04 




N160 






0.39±0.07 


239.1±0.1 


9.4±0.2 


3.81±0.08 


1.16±0.17 




C (»Pi ^ 




f 














N159W 






0.13±0.01 


237.6±0.1 


11.8±0.4 


1.59±0.05 


0.11±0.02 


6000 


N159E 






0.14±0.01 


235.5±0.1 


6.1it0.3 


0.87±0.03 


0.09±0.01 


7900 


N159S 






0.17±0.01 


238.1±0.1 


10.1±0.4 


1.90±0.06 


0.15±0.02 


1000 


N160 






0.14±0.01 


239.2±0.1 


2.4±0.1 


0.34±0.02 


0.10±0.02 


21700 


Rl 






0.06±0.01 


238.3±0.3 


9.1±0.7 


0.55±0.04 


0.14±0.02 


17800 


R2 






0.09±0.01 


236.6±0.2 


7.7±0.6 


0.73±0.05 


0.12±0.02 


1700 



''Errors are statistical and ±1(7. Systematic calibration uncertainty is 1<t ~ ±10%. 

''Values and errors derived from Gaussian fit. 

'^Ratio of transition intensity to that of ^■^CO(J = 1 — > 0), convolved to the same beam, in K km . 
To obtain ratios for intensities in erg s~^ cm~^ sr~^ multiply by (u/i/io)^. 

''Errors are ilu, and include 10% la calibration uncertainty in both lines. 

''[CIIJ/CO ( J = 1 — » 0) line intensity ratio in a 4' beam, in erg cm~^ sr"'^. 

* Measured in a HPBW=4' beam. 
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Table 4. FIR Parameters for the N159/N160 Complex 



Identifier Seo ^ Sioo ^ ^ T^^.^ ^ - Eq. Sp. Type ^ 

(Jy) (Jy) (K) (106 Lo) (ZAMS) 



N159W 


204.9 


671.8 


3.28 


30 


3.40 


5x05V 


N159E 


179.6 


393.7 


2.19 


34 


1.86 


3x05V 


N159S ^ 


6.6 


39.5 


6.01 


25 


< 0.24 


<06V 


N160 


393.1 


580.1 


1.48 


40 


2.74 


4x05V 


N160E 


48.5 


76.1 


1.57 


39 


0.36 


05.5V 



^Point source fiux density at 60 and 100 fxm derived from IRAS 
HIRES data. 

''Dust temperature derived from the ratio of Sioo to Seo, assuming 
the dust emits as a graybody with emissivity exponent (3=1. 

^FIR luminosity computed according to Lonsdale (1985), assuming 
a dust emissivity exponent P = 1 and D = 52 kpc. 

'^Equivalent spectral type and multiplicity (Panagia 1973), assuming 
all starlight is absorbed by dust and reemitted in the FIR. 

'^There is no infrared source associated with N159S, therefore the 
IRAS fluxes are upper limits. 
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Table 5. Radio Continuum Sources in the N159/N160 Complex 



Identifier 


R.A. 


Dec. 


Sl.42 ^ 




Sp. Type ^ 




(B1950) 


(B1950) 


(mJy) 


(pc cm~^) 


(ZAMS) 


N159-I 


5^'40™29!3 


-69°46'21" 


480 


143 


1.5x05V 


N159-II 


5''40""05!2 


-69°47'29" 


440 


139 


1.4x05V 


N159-III 


5h4Qm22S4 


-69°46'37" 


336 


127 


1.3x05V 


N159S 






< 30 


< 57 


<06V 


N160-I 


5'^40"07!8 


-69°40'17" 


802 


170 


1.7x05V 


N160-II 


5'^40°^47!3 


-69°41'42" 


169 


101 


05V 


N160-III 


5^40""35!1 


-69°41'49" 


126 


92 


05V 



^Point source flux density at 21 cm. 
''H II region excitation parameter. 

'^Equivalent spectral type and multiplicity, derived using U computed 
by Panagia (1973). 

^No radio continuum detected down to la ^ 10 mJy. 



